function [ xdot ] = model3(t, x)
%MODEL Summary of this function goes here
%   Detailed explanation goes here

global K productionRateFnc;

xdot = zeros(2, 1);
xdot(1) = productionRateFnc(t, x) - x(2);
xdot(2) = 6 * x(1) - K * productionRateFnc(t, x);

end

